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The distribution of cosmic-ray sources in the Galaxy, 7-rays and 
the gradient in the CO-to-H 2 relation 
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Abstract. We present a solution to the apparent discrepancy between the radial gradient in the diffuse Galactic 
7-ray emissivity and the distribution of supernova remnants, believed to be the sources of cosmic rays. Recent 
determinations of the pulsar distribution have made the discrepancy even more apparent. The problem is shown to 
be plausibly solved by a variation in the VKco-to-A (Eh) scaling factor. If this factor increases by a factor of 5-10 
from the inner to the outer Galaxy, as expected from the Galactic metallicity gradient, we show that the source 
distribution required to match the radial gradient of 7-rays can be reconciled with the distribution of supernova 
remnants as traced by current studies of pulsars. The resulting model fits the EGRET 7-ray profiles extremely 
well in longitude, and reproduces the mid-latitude inner Galaxy intensities better than previous models. 
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1. Introduction 

The puzzle of the Galactic 7-ray gradient goes back to 
the time of the C OS-B satellite (|Bloemen etail Il98fj 
IStrong et all Il988f) : using HI and CO surveys to trace 
the atomic and molecular gas, the Galactic distribution 
of emissivity per H atom is a measure of the cosmic-ray 
(CR) flux, for the gas-related bremsstrahlung and pion- 
decay components. However the gradient determined in 
this way is much smaller than expected if supernova rem- 
nants (SNR) are the sources of cosmic rays, as is gen- 
erally believed. This discrepancy was confirmed with the 
much more precise data from EGRET on the COMPTON 
Gamma Ray Observatory, even allowing for the fact that 
inverse-Compton emission (unrelate d to the gas) i s mor e 
important than originally supposed (Strong et all 120001. 
A possible explanation of the small gradient in terms of 
CR propagation, involving radial var iations of a Galactic 
wind, was recently put forward by iBreitschwerdt et alJ 
d2002t) . 

However the derivation of the Galactic distribution of 
SNR, commonly based on radio surveys, is subject to large 
observational selection effects, so that it can be argued 
that the discrepancy is not so serious. But other trac- 



ers of the distribution of SNR are available, in particular 
pulsars; the new sensitive Park es Multibeam s urvey with 
914 pulsars has been used by lLorimerl l)2004h to derive 
the Galactic distribution, and this confirms the concen- 
tration to the in ner Galaxy. Fig . ^ compares the pulsar 
distribution from lLorimerl l|2004jl with a CR source distri- 
butio n which fits the EGRET 7-ray data ijStrong et all 
2000). If the pulsar distribution indeed traces the SNR, 
then there is a serious di screpancy with 7-rays. The d is- 
tribution of SNR given bv lCase fc Bhattacharval l|l998h is 
not so peaked, but the number of known SNR is much 
less than the number of pulsars and the sys tematic ef- 
fects very difficult to account for l|Gree But even 
this flatter distribution is hard to reconcile with that re- 
quired for 7-rays. Another, quite independent, tracer of 
the SNR distribution is the 1809 keV line of 26 Al; whether 
this originates mainly in type II supernovae or mass- 
sive stars is not important in this context, since both 
trace star-formation/SNR. The COMPTEL 26 Al maps 
l|Knodlseder et a,llll999HPlfischke et a,ll l200lj) show that 
the emission is very concentrated to the inner radian of the 
Galaxy. T he density of free elect rons shows a similar dis- 
tribution l|Cordes fc Lazioll2003j) . The 26 Al measurements 
are not subject to the selection effects of other methods; 
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Fig. 1. CR source density as function of Galacto centric 
radius R. Dotted: as used 
line: based on pulsars ijLorimerl 



Strong et all l)2000ft . solid 
20041) as used in this work, 
vertic al bars: SNR data points from lCase fc Bhattacharval 
(1998). Distributions are normalized at R = 8.5 kpc. 



although they have their own uncertainties, they support 
the type of distribution which we adopt in this paper. 

A major uncertainty in the models of diffuse Galactic 
7-ray emission is the distribution of molecular hydro- 
gen, as traced by the integrated intensity of the J = 1— 
transition of 12 CO, Wgg. Gamma-ray analyses have 
in fact provided one of the standard values for the 
scaling factor 1 Xgo — N(H2)/Wco', with only the as- 
sumption that cosmic rays penetrate molecular clouds 
freely, the 7-ray values are free of the uncertainties of 
other methods (e.g. those based on the assumption of 
molecula r cloud virializati o n) . H o wever previous analy - 
ses, e.g. IStrong fc Mattel l)l996h . iHunter et aD (|l997|) . 
IStrong et alJ l)2000|h have usually assumed that Xgo is 
independent of Galactocentric radius R, since otherwise 
the model has too many free parameters. But there is now 
good reason to believe that Xgo increases with R, both 
from COBE/DIRBE studies <|Sodroski et all h 99,4 17997^ 
and from the measurement of a Galactic metallicity gradi- 
ent combined with the strong inver se dependence of X go 
on metallicity in external galaxies l|lsrael ll997L 12000^) . A 
rather rapid radial variation of Xgo is exj 



pected, based on 



a gradient in [O/H] of 0.04-0.07 dex/kpc faou et aHl200C ; 
Deharyeng et all l200(i iRolleston et all l2000l ISmartt , 



200 it lAndrievskv et all 12002ft and the dependence of 
Xgo o n metallicity in ext ernal galaxies: log Xgooc —2.5 
[O/H] lllsraell I1997L l2000lh giving AT co cx io(-° 14±0 04 )«, 
amounting to a factor 1.3-1.5 per kpc, or an order of 
magnitude between the inner and outer Galaxy. 2 A less 
rapid dependence, log Xggoc-1.0 [O/H], was found by 



units: molecules cm 2 / (K km s 

H 



2 The values given by llsraell Jl997l l200d) include the ef- 
fects of the radiation field, implicitly containing the radiation 
field/metallicity correlation of his galaxy sample. Xco is posi- 
tively and almost linearly correlated with radiation field, so the 
dependence of Xco f or constant radiation field is even larger: 
log Xcooc —4 [O/H] lllsraell 12000). By adopting the coefficient 
-2.5 we implicitly assume the same radiation/metallicity cor- 
relation within the Galaxy as over his galaxy sample. 



Fig. 2. Xgo as function of R. Dotted h o rizontal line, 
black: as used in IStrong k. Mattoxl l)l996|) : IStrong et all 
l|2000h : solid line, black: as used for -7-rays in this work; 
dashed, dark blue: from ISodroski et aD l)l995h : dash-dot, 
red: using m etallicity gra dient as described in the text, 
AT co oc Z- 2 b ftsraell EtOO). two lines for [O/H] = 0.04 and 
0.07 dex/kpc: dash- dot-dot, light blue: using Xgoct Z~ 1,n 
l|Boselli et all Eo02'l and [O/H] = 0.07 dex/kpc. The val- 
ues using metallicity are normalized approximately to 
those from the 7-ray analysis. 



iBoselli et alJ l|2002fl . w hich however still imp lies a signifi- 
cant Xgg>{R) variation. iBoissier et alJ l|2003|) also combine 
the metallicity gradient with Xgg.(Z) within individual 
galaxies, to obtain radial profiles of H2, and give argu - 
ments for the validity of this procedure. iDigel et alJ (^990) 
found that molecular clouds in the outer Galaxy (i?~12 
kpc) are undcrluminous in CO, with Xgo a factor 4±2 
times the inner Galaxy value. ISodroski et al 1 lll995Lll997l) 
derived a similar variation (logXco/10 20 = 0.1 2R- 0.34) 
when modelling dust emission for COBE data. IPak et all 
l|l998|) p redicted the physical o rigin f or a va riation of Xgp 
with Z. IPaoadopoulos et alJ l)2002() and IPaoadopoulosI 
ll2004l) discuss the physical state of this metal-poor gas 
phase in the outer parts of spiral galaxies (relatively 
warm and diffuse). Observations of H2 line emission f rom 
NGC 891 with ISO llValentiin fc van der Weri [l999) in- 
dicate a massive cool molecular component in the outer 
regions of this galaxy, supporting the trend found in our 
Galaxy. 

Fig. |2 illustrates some of the possible Xgo variations 
implied by these studies. For the cases where metallicity is 
used to estimate Xgo, the values are normalized approx- 
imately to the values used in the present 7-ray analysis, 
since we are only interested in comparing the variations of 
Xgo ■ From the viewpoint of 7-rays, the effect of a steeper 
CR source distribution is compensated by the increase of 
^CO Thus we might expect to resolve the apparent dis- 
crepancy in the source distribution, and improve our un- 
derstanding of the Galactic 7-ray emission. In this paper 
we investigate quantitatively this possibility. Note that the 
7-rays include major contributions from interactions with 
atomic hydrogen and from inverse Compton scattering, 
both of which are independent of Xgg] this means that 
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Fig. 3. Longitude profile of 7-rays for 1000-2000 MeV, 
averaged over \b\ < 5.5°. Vertical bars: EGRET data; lines 
are model components convolved with the EGRET point- 
spread function: green: inverse Compton emission, red: ir°- 
decay, light blue: bremsstrahlung, dark blue: total. 

the Xco variation has to be quite large to have a signifi- 
cant effect. 



2. Data 

The EGR ET and COMPTEL data are the same as de- 
scribed in lStrong et"afl l|2000L l2004a|) . The EGRET data 
consist of the standard product counts and exposure for 
30 MeV - 10 GeV, augmented with data for 10-120 GeV. 
The 7-ray point sources in the 3EG catalog ue have been 
removed as described in IStrong et ail (120001). The HI and 
CO data are as desc ribed in |Mosk^denkoetal and 
IStrong et alJ l)2004al) : they consist of combined surveys di- 
vided into 8 Galactocentric rings on the basis of kinematic 



ing models with data are given in IStrong et all ( 2004 
which the reader is referred. 



3. Model and Method 



to 



We u se the GALPROP program l|Strong et all [2000, 
2004a) to compute the models. GALPROP was extended 
to allow a variable Xqo{R) to be input. The distribu- 
tion of CR sources i s assumed to f ollow that of pulsars 
in the form given bv iLorimerl l|2004l ). as shown in Fig.^ 
The other parameters, in particular the CR nucleon and 
electron injection spectral shape and propag ation parame- 
ters, ar e taken from the "optimized model" of lStrong et alJ 
d2004al) . As before the halo height is taken as z h = 4 
kpc, and the maximum radi us R = 20 kpc . The is otropic 
background is as derived in IStrong et alJ (|2004tJ ). Since 
in this work we simply wish to demonstrate the possi- 
bility to obtain a plausible solution, we adopt a heuristic 
approach, adjusting Xco{R) to obtain a satisfactory solu- 
tion as shown in Fig. [2 The elec tron flux has been sca led 
down by a factor 0.7 relative to IStrong et al ] ll2004al) to 
obtain an optimal fit. 
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Fig. 4. Latitude profile of 7-rays for 1000-2000 MeV, av- 
eraged over 330° < / < 30°. Data and curves as in Fig.|3J 
The extragalactic background is shown as a black horizon- 
tal line. 

4. Results 

Figs.[3]andQ]show the longitude and latitude distributions 
for 1-2 GeV, compared to EGRET data. A rather rapid 
variation of Xqo is required to compensate the CR source 
gradient, but it is fully compatible with the expected vari- 
ation based on metallicity gradients and the COBE result, 
as described in the Introduction. The longitude and lati- 
tude fits are good except in the outer Galaxy where the 
prediction is rather low. One possible reason for this is 
that the CR source density does not fall off so fast beyond 
the Solar circle as given by the adopted pulsar distribu- 
tion, which has an exponential decay. Another possibility 
could be even larger amounts of H 2 in the outer Galaxy 
than we have assumed (see discussion in Introduction). 
We have chosen the range 1-2 GeV for the profiles since 
this is where the gas contribution and hence the effect of 
Xco is maximal. An exhaustive comparison of profiles in 
all energy ranges is beyond the scope of this Letter, but 
in fact the agreement is good at all energies. The larger 
CR gradient in this model has another consequence: the 
predicted inverse-Compton emission in the inner Galaxy is 
more intense at intermediate latitudes where the interstel- 
lar radiation field is stil l high; this is precisely the region 
where previo us models ijHunter et all Il997t IStrong et all 
120001120044 have had problems to reproduce the EGRET 
data. Fig.[S]shows the model spectrum of the inner Galaxy 
compar ed with EGRET data ; the fit is similar to that of 
models l|Strong et alll2004a|) with ad hoc source gradient 
and constant Xqo . The prediction is rather high above 20 
GeV, however the EGRET data are least certain in this 
range l)Strong et alll2004a|) . 

5. Discussion 

We have shown that a good fit to the EGRET data is 
obtained with the particular combination of parameters 
chosen. We can however ask whether the pulsar source dis- 
tribution combined with a constant A^co could also give 
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Fig. 5. Spectrum of inner Galaxy, 330° < I < 30°, |&j < 
5.5°. Vertical bars: EGRET data (red), COMPTEL data 
(green). Curves: predicted intensities; inverse Compton 
(green), 7r°-decay (red), bremsstrahlung (light blue), ex- 
tragalactic background (black), total (dark blue). 

a good fit if we reduce the CR electron intensity, to su- 
press the inner Galaxy peak from inverse Compton emis- 
sion. This can indeed reproduce the longitude profile in the 
inner Galaxy, but fails badly to account for the latitude 
distribution, since it has a large deficit at intermediate lat- 
itudes. Some variation of Xqo is therefore required. The 
suggested variation of Xco would have significant impact 
on the Galactic H2 mass and distribution. Warm molec- 
ular hydrogen in the outer parts of spiral galaxies that 
is not traced by CO emission may be detectable by the 
Spitzer observatory in 28 /im vibrational emission. These 
issues will be addressed in future work. 

6. Conclusions 

Two a priori motivated developments allow us to obtain 
a more physically plausible model for Galactic 7-rays, si- 
multaneously allowing a CR source distribution similar to 
SNR as traced by pulsars and an expected variation in the 
W / co-to--/V(H 2 ) conversion factor. Obviously the uncer- 
tainty in both the source distribution and Xqq are large 
so our solution is far from unique, but it demonstrates the 
possibility to obtain a physically-motivated model without 
resorting to an ad hoc source distribution. This result sup- 
ports the SNR origin of CR. The resulting model also gives 
improved predictions for 7-rays in the inner Galaxy at 
mid-latitudes. We have therefore achieved a step towards 
a better understanding of the diffuse Galactic 7-ray emis- 
sion. This result is important input to the development of 
models for the upcoming GLAST mission. This Letter is 
intended only to point out the potential importance of the 
effect. The next step will be a more quantitative analysis 
to derive Xqo(R) from the 7-ray data themselves. 
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